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We study the time evolution of a brane construction that is holographically dual to a strongly 
coupled gauge theory that dynamically breaks a global symmetry through the generation of an 
effective composite Higgs vev. The D3/D7 system with a background magnetic field or non-trivial 
gauge coupling (dilaton) profile displays the symmetry breaking. We study motion of the D7 brane 
in the background of the D3 branes. For small field inflation in the field theory the effective Higgs 
vev rolls from zero to the true vacuum value. We study what phenomenological dilaton profile 
generates the slow rolling needed, hence learning how the strongly coupled gauge theory's coupling 
must run. We note that evolution of our configuration in the holographic direction, representing the 
phyiscs of the strong interactions, can provide additional slowing of the roll time. Inflation seems 
to be favoured if the coupling changes by only a small amount or very gently. We speculate on how 
such a scenario could be realized away from N=4 gauge theory, for example, in a walking gauge 
theory. 



I. INTRODUCTION 

In recent years the AdS/CFT Correspondence [1-4] 
has emerged as a very powerful tool for studying strong 
coupling problems. It provides a weakly coupled gravita- 
tional/string description of strongly coupled gauge dy- 
namics allowing computation of many aspects of the 
dynamics. Originally the duality was proposed for the 
N = 4 super Yang Mills theory but it has been widely 
expanded to other theories over the past ten years. 
Amongst those theories for which gravity duals exist are 
several examples of strongly coupled gauge theories with 
fundamental quark fields and which break chiral symme- 
tries by the formation of a quark condensate as QCD 
does [5-8]. In this paper we wish to use these tools to 
make a first study of time dependence in such symme- 
try breaking theories. We expect strongly coupled gauge 
theories to undergo phase transitions between a high tem- 
perature chirally symmetric phase and a low temperature 
chiral symmetry breaking phase. Holography should for 
the first time allow us to analyse the evolution of such 
theories between these different phases. The simplest ex- 
ample of where such a transition might play an important 
role in the evolution of our Universe is in inflation and in 
this paper we will concentrate on that possibility. 

Inflation (see for example [9]) is now a key part of the 
standard cosmological model of the Universe, supported 
by many pieces of astrophysical data. The usual descrip- 
tion involves one or more scalar fields slow-rolling from 
an unstable point in a potential to the true vacuum. The 
origin for such a scalar field remains unclear though. In- 
deed, until evidence for supersymmetry is found in na- 
ture, fundamental scalar fields are formally unnatural in 
field theory as a result of the hierarchy problem (radia- 
tive corrections cause their mass to naturally grow to the 
Planck scale). This need be no obstacle to the paradigm 
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though since many strongly coupled systems have scalar 
order parameters describing their dynamics that could 
play the same role. For example, in QCD the gauge in- 
teractions dynamically generate a quark bilinear conden- 
sate (((?<?)) that breaks chiral symmetries - one may think 
of that bilinear as the expectation value of a composite 
scalar field. The main reason to work with fundamental 
scalars in cosmology is simply that we have had no tools 
to study strong coupling problems such as the condensa- 
tion in QCD. The AdS/CFT Correspondence for the first 
time allows us to analyse just such a scenario though. 

We will concentrate on a particular duality which we 
believe to be the simplest example of holography with 
fundamental quark fields [10-14]. We wish to stress that 
we do not consider the specific degrees of freedom of the 
theory too crucial - it is some strongly coupled gauge 
theory that generates quark condensates. We hope, in 
the spirit of AdS/QCD modcls[15, 16], that it reflects 
broad aspects of many strongly coupled systems. The 
specific gauge theory is constructed from the D3/D7 sys- 
tem 1 in type IIB string theory which we will describe in 
detail below. The theory is the large N c = 4 U(N C ) 
gauge theory with a small number of quark hypermul- 
tiplets. We will work in the quenched approximation 
[12] (appropriate when Nf -C N c ) which on the gravity 
dual side corresponds to treating the D7 branes as probes 
in the AdS metric generated by the D3 branes. There is 
a U(l) symmetry (a remnant of the SU(A) R-symmctry 
of the 3\f = 4 theory) which is broken when a quark con- 



Thc D3/D7 system has been used to construct an inflation model 
in [17] and subsequent literature. The motivation in those models 
is rather different since they are not holographic descriptions of 
a strongly coupled gauge theory but instead assumed to describe 
weakly coupled fields with the extra dimensions compactificd. 
The dilaton profiles and brane motions we will consider are very 
different being inspired by the gauge duality and the lessons we 
will seek to extract are for strongly coupled gauge dynamics. 
Another interesting AdS/CFT approach for a strongly coupled 
inflaton was reported in [18] 
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dcnsate forms [6] . Several mechanisms for triggering this 
condensation have been explored. The cleanest is when 
a background magnetic field is introduced [19, 20]. Run- 
ning of the coupling also causes quark condensation as 
has been shown in back-reacted dilaton flow geometries 
[6, 7]and models with a phenomcnologically imposed dila- 
ton profile [8]. The quark condensate can be determined 
in these models and an effective IR quark mass is gener- 
ated. The theories display a massless pion-like Goldstonc 
field and a massive sigma field (since we are at large N c 
it is stable) that is the effective Higgs particle. 

We will seek to learn about a strongly coupled counter 
part to a "new inflation" or "small field inflation" model 
[21]. One assumes a high temperature phase where the 
effective Higgs scalar's vev is zero and a low tempera- 
ture phase where the vev is non-zero. A second order 
phase transition should connect these phases so that the 
vacuum is left at the unstable zero vev point of the low 
temperature theory and then rolls to the true vacuum. 
For inflation the potential needs to be very flat near the 
origin so the roll takes a long time. Of course, other 
possible inflationary scenarios exist with multiple scalar 
fields or where the scalar expectation value is initially 
large relative to the vacuum value and it would be in- 
teresting to investigate these ideas holographically too in 
the future. 

In principle one could imagine taking a holographic 
model and tracking it through the finite temperature 
transition. On the gravity side temperature is introduced 
through a black hole in the AdS geometry. The D3/D7 
model we investigate has been shown to display both 
first order (eg with a magnetic field [19, 20]) and sec- 
ond order (eg with a magnetic field and chemical poten- 
tial [20]) symmetry restoration transitions. Ideally one 
would work in a time dependent background describing 
a shrinking black hole. In fact such geometries are known 
[22] and a first study of D7 branes moving slowly near a 
potential minimum in those geometries can be found in 
[23] . Dealing with such branes when they touch the black 
hole is hard though - we hope to study this problem in a 
future publication [24]. 

Our goal in this paper is more limited though and sim- 
ilar to the usual simplest field theory analysis of inflation 
models. We will assume that below the phase transition 
the theory is well described by the T = theory. We will 
invoke initial conditions that place the vacuum in the 
symmetric phase and then watch it roll to the true vac- 
uum that breaks the symmetry. These initial conditions 
assume the existence of a second order transition. 2 Our 
first results are numerical simulations of this roll in the 
theory where a magnetic field is inducing the symmetry 
breaking. As we said above this theory actually has a 
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first order thermal transition rather than a second order 
one - the scenario serves to demonstrate the formalism 
though. The results would be relevant to the phase tran- 
sition within a highly supercooled bubble. Here there is 
no fine tuning to make the potential particularly flat so 
the model would make a poor inflation model. It serves 
as a proof of the numerical techniques we use to study 
the problem and it is impressive that a relatively sim- 
ple computation can track such a transition in a strongly 
coupled gauge theory. 

One does not imagine that a magnetic field induced 
symmetry breaking is likely in a realistic model of in- 
flation. We therefore turn to a more phenomenological 
model in which we embed the D7 brane in an AdS geome- 
try but with a hand chosen and non-back reacted dilaton 
profile. Whilst this is not a completely kosher string dual 
we hope that it allows us to understand broadly the de- 
pendence of a strongly coupled theory's dependence on 
its running coupling. We used a similar approximation 
in [8] to study walking gauge dynamics and found agree- 
ment between the gravity dual description and the usual 
expectations for the coupling dependence of the chiral 
symmetry breaking. The model is sufficiently simple that 
we can understand how the running of the gauge coupling 
affects the potential for the effective Higgs mode. We at- 
tempt to engineer models with a slow roll between the 
symmetric and broken phases. We also show that part 
of our success in slowing the roll time is due to addi- 
tional dynamics in the holographic radial coordinate of 
the dual description. This dynamics reflects the strong 
coupling dynamics of the gauge theory which seems under 
some circumstacnes to favour inflation. From this anal- 
ysis we conclude that gauge theories displaying a small 
increase in the gauge coupling, preferably over a wide 
energy regime, are liable to make good inflation models. 

Finally we briefly speculate on how IR conformal fixed 
points may well be common in strongly coupled gauge 
dynamics and suggest a wider set of theories that might 
display running of the coupling like that we find gives 
inflation. Obvious examples are walking gauge theories. 

II. INFLATION 

Let us very briefly review the standard inflation model 
[9]. Based on the cosmological principle we assume the 
Friedmann-Robertson- Walker metric of the form 

/ dr 2 \ 

ds 2 = -dt 2 + alt) 2 + r 2 (d9 2 + sin 2 9d(j?) , 

\ 1 — kr z J 

(1) 

and the matter energy momentum tensor of a perfect 
fluid Tjf = diag(— e,p,p,p). For simplicity we will work 
in the k = case. Then the general Einstein equations 
describing the expansion of the Universe are reduced to 
two equations: 

H 2 = — e : Friedmann equation , (2) 
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£ + 3H(e +p)=0 : Fluid equation , (3) 

where H(t) = a/a is the Hubble parameter and the scale 
factor a(t) is readily determined from H(t) as 



a(t) = o(0) exp ( y ff(t')(2t' 



(4) 



Equations (2) and (3) can be combined to form a so 
called "acceleration equation" 



o 1 / . 

= -77 £ + 3p , 
a o 



(5) 



from which the condition for inflation can be expressed 
in terms of the energy momentum tensor. 



Inflation <^ d> <^> P < ~o 



(6) 



Thus a natural question for inflation is asking what kind 
of matter and dynamics can generate sufficient negative 
pressure (6). 

For example let us consider a scalar field. 



-9^ + Sm 



where 31 is Ricci scalar and 



Sm — 



(7) 



(8) 



in unit of 8nG = 8nmp =1. nip is the Plank mass. 
From a scalar matter action the energy-momentum ten- 
sor of a homogeneous field ip(t) can be obtained as 

e = -T ° = \f + V(tp), v = Tt = \p % - Vfr) . (9) 
The inflation condition (6) is rephrased as 



P< 



<^> tp 2 < V(if) . 



(10) 



Whenever the kinetic term is small compared to the po- 
tential energy there will be inflation. Then by analogy 
to mechanics we may say tp is rolling slowly. 

Let us apply this slow roll condition to the equations 
of motion (2) and (3). 



more if tp is sufficiently smaller than the other terms in 
(12) then V(<p) can be considered as constant, say Vq, for 
a long time. Thus the scale factor (4) in the slow rolling 
range reads 



a(t) ~ a(0)e m = o(0) 



(13) 



where we have reinserted the Planck mass. This is the 
de Sitter limit zeroth order approximation of an inflating 
Universe, which we will adopt in this paper. To have a 
graceful exit from inflation the Hubble parameter must 
be time dependent and fall to zero at late times. In an 
inflating regime the space is equivalent to dc Sitter space 
with a cosmological constant A = 8irrrip 2 Vo, where Vb is 
interpreted as the vacuum energy density. The amount 
of inflation is specified by the number of c-folds given by 



N, 



, a ( f e) 
= l0S ^) 



3 TOp 



(14) 



where t e is the time when inflation ends (when slow roll 
conditions are violated) starting from t = 0. One needs 
N e to be more than 60 phenomenologically. 

Note that a larger H yields a more inflationary evo- 
lution. H is nothing but a friction coefficient when we 
interpret (12) as an equation for the classical particle tra- 
jectory {<p{t)) under the potential V. Thus it is natural 
that a larger friction induces a slower rolling of the par- 
ticle. 

In this paper we will replace the scalar field sector with 
a strongly coupled gauge theory. The gauge dynamics 
will generate a wine bottle shaped potential with a non- 
zero condensate of a quark bilinear at the minimum. We 
will phenomenologically adjust the running of the gauge 
theory's coupling constant to control the shape of the po- 
tential. We hope to learn how a gauge theory's coupling 
must run to generate inflation. Clearly our task is to 
compute the stress energy tensor (e and p) of a strongly 
coupled gauge theory so we can substitute them into the 
Fricdmann and fluid equations above. We must find the 
equivalent of the Euler Lagrange equation for the scalar 
field tp in the inflating background for the time evolution 
of the strongly coupled gauge dynamics. To do this we 
will turn to Gauge/Gravity duality. 



III. HOLOGRAPHIC DESCRIPTION OF A 
STRONGLY COUPLED GAUGE THEORY 



H 



1 



1 /l 



3 V 2 



-tp 2 + V(tp) 



dV 

tp + 3Htp + — = 
dp 



(11) 
(12) 



where the fluid equation (12) is nothing but the equation 
of motion for tp, which can be obtained from (8). 

If the scalar field stays in some part of the potential 
with tp 2 < V then H - V(tp(i))/3 from (11). Further- 



A. Brane Construction 

First let us review the gravity dual description of the 
symmetry breaking behaviour of our strongly coupled 
gauge theory [5-8]. To begin with we will work in a time 
independent flat space (i.e. there is no inflation here and 
a(t) = 1). 

Dp-branes are p dimensional membrane- like objects on 
which the ends of open strings can be fixed. The weak 
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FIG. 1: A sketch of the D3/D7 construction that our model 
is based on. 3-3 strings generate the gauge theory whilst 3-7 
strings are quark fields. 



coupling picture for our D3/D7 set up is shown in Fig 1 - 
there are iV D3 branes and the lightest string states with 
both ends on that stack generate the adjoint representa- 
tion fields of the N = 4 gauge theory. Strings stretched 
between the D3 and the D7 are the quark fields lying 
in the fundamental representation of the SU(N) group 
(they have just one end on the D3) . 

In the strong coupling limit the D3 branes in this pic- 
ture are replaced by the geometry that they induce [1-4]. 
We will consider a gauge theory with a holographic dual 
described by the Einstein frame geometry AdS$ x S 5 



ds 2 



1 

guv 



r 2 R 2 

-jpdxl + — {dp 2 + p 2 dtt\ + dwj + dwl) 



(15) 

where we have split the coordinates into the 2:4(2:3+1) of 
the gauge theory, the p and il 3 which will be on the D7 
branc world- volume and two directions transverse to the 



D7, W5,u>6- The radial coordinate, 



= P 2 + w\ +w 2 , 
corresponds to the energy scale of the gauge theory. The 
radius of curvature is given by R 4 = Airg^yNcOt' 2 with 
N c the number of colours. In addition we will allow an 
arbitrary dilaton as r — > to represent the running of 
the gauge theory coupling 



(16) 



where the function j3(r) — > 1 as r — > 00. The r — > 00 
limit of this theory is dual to the N = 4 super Yang-Mills 
theory and g^y is the constant large r asymptotic value 
of the gauge coupling. 

We will introduce a single D7 brane probe [12] into 
the geometry to include quarks - by treating the D7 as 
a probe we are working in a quenched approximation al- 
though we can reintroduce some aspects of quark loops 
through the running coupling's form if we wish (or know 
how). This system has a U(l) symmetry acting on the 
quarks, corresponding to rotations in the w^-Wq plane, 
which will be broken by the formation of a quark con- 
densate. 

In the true vacuum at zero temperature (T = 0) the 
brane will be static. We must find the D7 embedding 



function e.g.w^(p) = L(p) with wq 
in Einstein frame is given by 



0. The DBI action 



S D 7 = -T 7 fd*Ze*y/-detP[G\ ab 



(17) 



-T 7 J d 4 x dp p 3 /Vl + L' 2 



where P[G] a b is the pull back of the background metric 
onto the D7 and L' = d p L(p). T 7 = 1/(2tt)V 4 and % = 
2w 2 T 7 /g 2 JV when we have integrated over the 3-sphere 
on the D7. The equation of motion for the embedding 
function is therefore 



Vl + L' 2 _ 



2LpVl + £' 2 



d/3(r 2 
dr 2 



(18) 



The UV asymptotic of this equation, provided the dilaton 
returns to a constant so the UV dual is the N = 4 super 
Yang-Mills theory, has solutions of the form 



L = m 



(19) 



where we can interpret m as the quark mass (m q = 
m/2-Ko 1 ) and c is proportional to the quark condensate. 

The embedding equation (18) clearly has regular so- 
lutions L = m when g\ M is independent of r - the flat 
embeddings of the N = 2 theory [12, 13]. Equally clearly 
if d(3(r 2 )/dr 2 is none trivial in L then the second term 
in (18) will not vanish for a flat embedding. 

There is always a solution L = which corresponds 
to a massless quark with zero quark condensate (c = 0). 
In the high T phase this is the true vacuum. In the 
symmetry breaking low T geometry this configuration is 
a local maximum in the potential. 

Note that in the particular case when 




(20) 



the DBI action for the D7 brane is that of the D3/D7 
system with a background magnetic field B. This model 
has been extensively studied in [19, 20]. The action man- 
ifestly grows as one approaches L = p = so the D7 
brane is repelled from that point. In Fig 2 we plot the 
D7 embedding in the magnetic field case. 

An interesting phenomcnological case is to consider a 
gauge coupling running with a step of the form 



/3 = J 4 + 1-Atanh[r(r-A)] 



(21) 



Of course in this case the geometry is not back reacted 
to the dilaton and the model is a phenomcnological one 
in the spirit of AdS/QCD. This form introduces confor- 
mal symmetry breaking at the scale A = A/27ra' which 
triggers chiral symmetry breaking. The parameter A de- 
termines the increase in the coupling across the step. If 
the coupling is larger near the origin then again the D7 
brane will be repelled from the origin. The parameter T 
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FIG. 2: The coupling flow (16) (top) and the D7 brane em- 
beddings/quark self energy (bottom). 



spreads the increase in the coupling over a region in r of 
order T _1 in size - the effect of widening the step is to 
enhance the large p tail of the D7 embedding. 

We again display the embeddings for some particular 
cases in Fig 2. Note that we have chosen parameters here 
that make the potential difference between the symmetric 
and symmetry broken phases the same in each case. This 
is crucial to ensure that we are comparing models that 
will generate the same cosmological constant and hence 
the same rate of inflation when the quark condensate is 
zero. The vacuum energy is given by the DBI action 
evaluated on the solution. In fact this energy is formally 
divergent corresponding to the usual cosmological con- 
stant problem in field theory. As usual we subtract the 
UV component of the energy - we do this by subtract- 
ing the energy of the symmetry breaking, lowest energy 
embedding in each case. 

The symmetry breaking of these solutions is manifest. 
The U(l) symmetry corresponds to rotations of the so- 
lution in the w^-wq plane. An embedding along the axis 
corresponds to a massless quark with the symmetry un- 




FIG. 3: The approximate potential, V = Sm/Tj, for the 
quark condensate c, from the holographic model for the run- 
ning couplings and D7 embeddings of Fig 2. 



broken (this is the configuration that is preferred at high 
temperature and it has zero condensate c). The vacuum, 
curved configurations though map onto that case at large 
p (the UV of the theory) but bend off axis breaking the 
symmetry in the IR. 

One can interpret the D7 embedding function as the 
dynamical self energy of the quark, similar to that emerg- 
ing from a gap equation [27]. The separation of the D7 
from the p axis is the mass at some particular energy 
scale given by p - in the 3\f = 2 theory where the embed- 
ding is flat the mass is not renormalized, whilst with the 
magnetic field or running coupling an IR mass forms. 



B. Approximate Potentials 

It is natural to want to plot the potential for the quark 
condensate using the holographic description. However, 
this is somewhat ambiguous. The embedding equation 
determines the D7 embeddings that correspond to the 
turning points of such a potential. In between these 
points one would need to find the minimum action con- 
figuration that falls off in the UV as c/p 2 for arbitrary c 
and has L' = in the IR. This can't be done by a simple 
numerical shooting because such configurations are not 
solutions of the Euler Lagrange equations. A reasonable 
way to get qualitative results though is as follows - we 
simply take the vacuum embedding solution and plot the 
potential as a function of that embedding multiplied by 
an arbitrary constant. This gives embeddings for all L(0) 
and all values of c. The case when we multiply by zero is 
of course the central maximum of the wine bottle poten- 
tial. We will assume that this form is appropriate for all 
values of L(0) or c. This then lets us plot the potential for 
each of our models above - the potential automatically 
takes the correct value at the maximum and minimum. 
See Fig 3. 

Although this procedure is somewhat adhoc it allows 
a first rough understanding of the potential shape. We 
can see that the potential for the case of a magnetic field 
induced symmetry breaking is rather steep around the 
origin. The step function ansatz for j3 is gentler there 
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and by decreasing A the curvature is further reduced. 
Note that as one decreases A one must increase A in or- 
der to keep the value of the potential when c = equal. 
In the extreme low A limit the model is characterized by 
that potential value being much less than the character- 
istic scale of the running A. The quark condensate also 
grows in this limit. Comparing the blue curve and the 
red curve shows that varying T also leads to a potential 
with a flatter and longer period near the origin. This is 
encouraging since for an inflation model we would want 
the potential to be as flat and as extended as possible 
around the origin. 

The lowest A curve (Orange) shows a new minimum 
at zero and an intermediate maximum in the potential 
but this is an artefact of the crude approximation be- 
ing made here. Were the maximum to really form then 
there would be a new solution to the embedding Euler- 
Lagrange equation for the D7 brane - we have checked 
and no such solution exists. We are simply describing 
the non-vacuum configurations incorrectly. Fig 3 is we 
stress just to show qualitatively that a Higgs-like poten- 
tial exists in the theory. Indeed we will see below that 
the time evolution of a rolling configuration depends cru- 
cially on the behaviour in the holographic directions of 
the description of the gauge theory which these plots are 
blind to. To do a better job of understanding the dynam- 
ics of these models near the origin we turn to studying 
the full time dependent problem of a moving D7 brane. 



IV. HOLOGRAPHY OF INFLATION 

A. Time Dependent Holography 

We will now move on to look at time dependent evo- 
lution in the gauge theory described above. Strictly the 
AdS / CFT Correspondence describes our gauge theory in 
rigid flat space. To include inflation we will instead as- 
sume a background holographic geometry given by 



ds 2 



1 

3uv [ 



R 2 



(gttdt 2 + gijdx l dx-') 



(dp 2 + P 2 dn 



dw 2 b 



dw 2 ) 



(22) 



where gu = —1 and gij = a(t) 2 Sij. For this to be valid 
even as an approximation we must formally assume that 
a(t) is growing much slower than the typical time for 
equilibration in the N = 4 Yang-Mills plasma. How then 
can we hope to study inflation with exponential space- 
time growth? The formal answer is that we will study 
inflation driven by the flavour sector of our gauge theory. 
The contribution of the D7 branes to the full vacuum en- 
ergy of the field theory is 1/N C suppressed relative to the 
glue dynamics contributions. We will be able to study 
the theory where this perturbation to the vacuum en- 
ergy drives inflation but formally in a limit where the 



inflation is still slow relative to the 3\f = 4 equilibration 
time. Our inflation is that relevant for the quark fields. 
This provides a formal justification for what follows. 

In fact for a true model of inflation one would want 
exponential expansion relative to all time scales in the 
problem. The reader can therefore choose to view (22) 
as a phenomenological model. In practice for the results 
we wish to extract - what dilaton profile will generate 
slow rolling behaviour - one could just work in static X4 
flat space. However, it is useful to have the damping 
term in the scalar equation of motion in (12) to remove 
late time oscillations - this provides another justification 
for using (22). 

The DBI action in such a geometry is given by 



''DBI 



(23) 



where 



L = -T7R 4 



dpPp* 



fJ 



00^2 



(p 2 + L 2 y 



1 + L' 2 (24) 



Here g M „ are boundary field theory metric components 
and L = u>5(t, p) with ojq = 0. L' = d p L and L = d t L. All 
variables (p, L, t, V, A) in the integral (24) are rescaled by 
R and dimensionlcss. We will reinstate R when needed. 
Note that the 4D effective action derived from the DBI 
action can be written in a covariant form with the correct 
measure. 

The integral form of the 4d Lagrangian (24) may be 
understood as the effective Lagrangian of some field the- 
ory quantity after integration over the extra direction 
p. Since the bulk embedding dynamics is closely related 
to the spontaneous symmetry breaking in the boundary 
field theory, it is natural to associate an effective degree 
of freedom to an order parameter, the condensate c(t). 
Thus we may consider (24) as an inflaton model of a com- 
posite scalar field c(t) with some potential V(c). The 
quark condensate c can be extracted numerically from 
the asymptotic large p form of the solution (19). Even 
though it's not straightforward to read the form of V(c) 
from (24) it is not an obstacle to the study of inflation, 
since s and p can be computed from (24) and we can ap- 
ply the inflation condition (6). Furthermore our potential 
for the condensate has its origin in (24), so is determined 
in principle and not put in by hand. 

We can obtain e and p by computing the expectation 
value of the stress energy tensor of the 3+ Id field theory 
from 



2 SS 



M 



-g Sg^ v 



(25) 



According to the AdS/CFT Correspondence the gravita- 
tional action is the Masterfield of the field theory [1-4] 
with the boundary values of gravitational fields playing 
the roles of the sources. Here we have explicitly shown 
the behaviour of the gravity action on the boundary met- 
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ric so we may simply compute. We find 



dimcnsionlcss. 



e{t) = T 7 R 4 / dpp 3 /3 



l + L 



1 + L' 2 



p(t) = -T 7 i? 4 / dpp 6 p 




L 2 

5W 



- £o 



(26) 



L' 2 - 



I? 



(p 2 + 



-Po (27) 



where we have renormalizcd by subtracting Eq,pq, which 
are the values of the integrals with the asymptotic static 
symmetry breaking solution for L, say L s . Examples of 
L s are shown in Fig 2. There is no explicit a(t) depen- 
dence in e and p, but its information is encoded in L(t, p) 
since L(t,p) is a solution of equation of motion of (23) 
which includes a(t) in \J —g. 

Note that p(0) = -e(0) if L = 0. So if an initial 
condition has L ~ then the Universe will start inflat- 
ing (p < — e/3) regardless of an initial L'. However we 
know that eventually the embedding will asymptote to 
L s because of the wine bottle shaped potential, which 
implies e — s- and p — > 0. Thus we expect e is decreas- 
ing and p is increasing with time. If p increases faster 
than e then there will be a time when e + Zp = 0, ending 
the inflation. How long it takes depends on the initial 
configuration L(0,p) and the parameters of /3, the dila- 
ton profile. According to small field inflation models it 
is natural to start with L = 0, which is a local maxi- 
mum configuration and the symmetric minimum at high 
temperature. Then the inflation time or the number of 
e- folds (iV e ) can be studied MS M function of coupling (3. 
In our setup inflation always happens to some degree and 
the issue is the amount of inflation. 

Since we know e and p we can proceed with the Fried- 
mann (2) and fluid (3) equations. Let's start with the 
Friedman equation. 



" 2 = r- 



(28) 



where e is (26). This is very complicated, e has H(t) in 
a(t) which is enters in the equation of motion for L. Also 
the Fricdmann equation is coupled to the fluid equation. 

To make the computation tractable we start with the 
zeroth order slow rolling approximation by assuming e is 
almost constant for a long enough time, namely e = Vo- 
lt gives us a simple solution for a(t) 



a(t) 



(29) 



where we set a(0) = 1. This solution must be consistent 
with the fluid equation, which is equivalent to the equa- 
tion of motion for L. Thus we will solve the equation of 
motion resulting from (23) with a constant Vq (29) and 
then plug the solution back into (26). If the calculated 
e changes slowly enough (slow rolling), then our solution 
is self-consistent. In our numerical computation e and p 
are always rescaled as e — > eT^R 4 and p —> pT^R 4 so are 



B. Rolling in a B field 

As a first example of our formalism we will consider 
the case of magnetic field induced symmetry breaking. 
We do not expect such a model to be well suited to a 
realistic inflation model because the typical curvature of 
the potential (Fig 3) for the condensate c is large. There 
is no particular fine tuning in the model. We will be 
able to track the time evolution of the brane configura- 
tion though from a symmetric to a symmetry breaking 
vacuum. 

We will first study the time dependence of the model 
in a constantly inflating Universe i. e. with fixed cosmo- 
logical constant H. In reality H should be determined 
by the D7 energy density through the evolution (H de- 
pends on the choices of N c and the 'tHooft coupling - we 
will show some generic numerical results). In the early 
inflating stages of the evolution we are most interested 
in, assuming a constant H is sound. For the evolution 
near the true vacuum H should go to zero and the brane 
would be expected to oscillate about the true vacuum 
configuration. 

Most of the solutions we will show will remain highly 
damped in the late time regime which allows us to nu- 
merically test that the configuration indeed ends on the 
static vacuum D7 embedding. We will though show some 
results for H = 0, where those oscillations are observable, 
shortly. The high damping is due to our unphysical con- 
stant H at all time. In a full computation it would be 
time dependent and vanish at late time. If we could solve 
with a time dependent H we would find the oscillation 
around the true vacuum at late time as hinted in the case 
H = 0. 

In particular we will start with the initial conditions 



L(p,t = 0) = 0, 



L(p,t = 0) = ve' p 



(30) 



where L = dtL. The initial speed and p dependence of 
this ansatz is not picked for any deep reason but is just 
illustrative of some initial condition that initiates the roll 
down to the potential minimum. We have checked that 
none of our results are qualitatively changed by varying 
for example the width of this initial condition. We will 
typically pick v to be very small so that the roll time from 
the peak of the potential is quite long. This also ensures 
that inflation happens at early times. The early time 
inflationary period in the plots below are the result of this 
fine tuned small initial condition and the large damping 
term and not a sign that the potential is particularly flat 
near the origin. 

A slow early roll will mean we can study the dynamics 
more easily numerically and understand better how to 
lengthen that early roll period. Through the roll we must 
ensure that L'(0,t) = and that L(oo,t) = 0. The 
evolution can in fact be followed with these boundary 
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L(l,p), with B=35.6 




FIG. 4: We show the time evolution of the D7 brane with 
H = y70/3, B = 35.6 and initial velocity parameter v — 
0.00001. The transition from the flat embedding (chirally 
symmetric phase) to the curved embedding (chirally broken 
phase)is apparent. 



\ . 

20 40 60 80 i 100 120 



' 

FIG. 5: e + 3p(dotted) and e(solid) plotted against time for 
the time evolution of the configuration shown in Fig 4. 



conditions using the in built numerical partial differential 
equation solver in Mathematics. 

In Fig 4 we show a sample plot of the numerical 

evolution. The figure is for an inflating Universe with 

H = ^/70/3, B = 35.6. We show a three dimensional 
plot for the time evolution of the D7 embedding - at early 
times the D7 is flat, L = 0, but as soon as the kinetic 
energy of the brane begins to grow as it experiences the 
curvature of the potential, it rapidly transitions to the 
vacuum embedding with symmetry breaking. 

Fig 5 is a plot of e + 3p and e versus time. e(t) sa —p{t) 
at all time, e shows a plateau until t ~ 80 and changes 
abruptly, while e + 3p is negative. Thus our slow roll 
approximation is valid up to t ~ 80. Even though the 
solution after t ~ 80 is beyond our approximation we read 
off the inflation ending time as the time when e + 3p = 
as an estimate, which is t e ~ 93. The end time will 
depends on B or H, which are related. 

We stress again that we have taken the Hubble param- 
eter H constant through the brane motion in the results 
just presented. With this unphysical late time damping 
it can be seen that the solutions precisely match on to 
the time independent vacuum configuration at large time. 



c 




FIG. 6: We show the time evolution of the quark conden- 
sate with H — 0, B = 30 and initial velocity parameter 
v = 0.00001. The D7 brane oscillates from the flat embedding 
(chirally symmetric phase) to the curved embedding (chirally 
broken phase) above and below the axis. 



To show the oscillatory behaviour one expects when the 
Hubble parameter is not present we can also solve for a 
similar configuration with H = 0. We show a plot in Fig 
6 of the time dependence of the quark condensate c in 
such a scenario - c is extracted from the large p depen- 
dence of the solution through (19). With no damping at 
all the D7 moves to approximately the vacuum config- 
uration, then overshoots, returns to the flat embedding 
before moving below the axis, etc. The oscillatory be- 
haviour is clear and can be followed through many cy- 
cles. Note that there are "fine wrinkles" near c = in 
this plot. These are due to a peculiarity of the magnetic 
field induced symmetry breaking - in particular there are 
an infinite set of meta-stable vacua near c = in this 
theory as explored in [28] . We see their influence on this 
motion although they play no particular role in our anal- 
ysis here. Our phenomcnological dilaton profiles below 
do not generate such structure. 

Let us briefly return to compare these simulations to 
what would be needed for inflation in our Universe. The 
energy density (23) is naively 

e = f 7 R 4 E(B), (31) 

where E(_B) is a number obtained by the numerical in- 
tegration in (24). Note the choice of the magnetic field 
here introduced the intrinsic scale of symmetry breaking. 

e is also a measure of the symmetry breaking scale of 
the theory and in flat space we could set it to one as 
the defining scale in the problem. When we include the 
inflationary dynamics through the damping term in (24) 
the Planck scale also enters. The ratio of the Planck 
scale to the fourth root of £ is a free parameter. We 
will associate the Planck mass with some distance in the 
holographic p direction, j p . The physical Planck mass in 
the coordinates with which we compute (rescaled by R 
as discussed under (24)) is thus 

m P = % . (32) 
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It is then sensible to express the energy density in Planck 
units so 



where 



E(B)- 



e = eo m p 



R s 



(33) 



, (34) 



(2tt)3 7 4 



with \ t fj = g^jyN . Note this is a' independent as it 
should be since R 8 /a' 4 = (4ngfj V N) 2 . 

The roll times t e we have computed are in dimension- 
less units and should be written in Planckian units too. 

t e m P = Rt e m p = -y p t e . (35) 

The scale factor in the metric is therefore given by 

/ 8ire f 



ait) 



Jit 



(36) 



Our simulations have been with J ^p7 P = -\/70/3. 

Note one can realize this for any value of eo - once so 
is fixed we can choose 7 p by setting the gauge theory 
parameters such as Xtn and N c . 

We conclude this section with the comment that it is 
quite remarkable that we can compute so straight for- 
wardly the time evolution of a strongly coupled gauge 
theory! 



C. Towards an Inflationary Dilaton Profile 

Our next goal is to look to extend the period of slow roll 
in the holographic models by phenomcnologically chang- 
ing the dilaton profile or gauge coupling's running. This 
is the equivalent of making the usual inflaton potential 
flatter around the origin. 

To examine this analytically let us first consider the 
roll in static space H = 0. We linearize the equation of 
motion about the initial symmetry preserving configura- 
tion L = 



L = 3p 3 L' 



At I. 

P L 



P 



I MO 



AT 



P -Q-L, 

Po 



(37) 



where /?o is the coupling at L = (for the step function 
form of the running we have /3q = A+l— Atanh[r(p— A)] 
and (3' = — Arsech[r(A— p)] 2 ). If we consider a truly flat 
configuration so L' = L =0 then only the last term on 
the right hand side contributes. Clearly the acceleration 
of the brane is localized around the point in p where 
there is a significant change in the coupling value. We 
show that this is indeed how the motion proceeds in Fig 
7 where we show a plot of the early time evolution of 
the brane in the presence of a step in the value of the 
dilaton. A bump grows around the point in p where the 
dilaton changes. We also plot the full evolution with a 



L(t,p) 




0.00015 

'o.oooio 

: 0.00005 
1 0.00000 



FIG. 7: The early time evolution of the D7 brane in the 
presence of a dilaton step with A — 3, A = 3. 24, V — 1 
and for the intial condition v = 0.00001. H = y/70/3. The 
acceleration occurs around the step. 



L(t,p) 




FIG. 8: The full time evolution of the D7 brane in the pres- 
ence of a dilaton step with A = 3, A = 3.24, V = 1 and for 
the initial condition v = 0.00001. H = ^70/3. 



cosmological damping parameter H = ij/70/3 showing 
the configuration move to the true vacuum in Fig 8. 

For the step configuration with large T, so we can con- 
sider it to be a sharp step, we have 



Hp) 



\ 3 ^jL(\)5(p-\) , 



(38) 



naively lowering A or T will reduce the acceleration. In 
fact though to keep the difference in the potential be- 
tween the symmetric and symmetry breaking configura- 
tions equal as we reduce A or T we must increase A. For 
the numerical values in Fig 3, increasing A wins so 7^4 A 3 
increases when A decreases. 

This analysis is overly naive though because as the D7 
brane evolves, as shown in Fig 7, L' and L will become 
important. Around the localized peak of the bump where 
L' = L is negative and the equation of motion will be 



Hp) 



A 



AT 

A + l' 



L(\) + \ 4 L"(\))5(p-\) . (39) 
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FIG. 9: Plots of e vs t for the three configurations [A — 3, 
r = 1, A = 3.24], [A = 1.8, r = 1, A = 4.325] and [A = 1.8, 
r = 0.5, A = 5.882]. v = 0.00001. H = ^70/3. 



FIG. 10: The condensate c vs t, for the three configurations 
[A = 3, T = 1, A = 3.24], [A = 1.8, r = 1, A = 4.325 ] and 
[A = 1.8, T = 0.5, A = 5.882]. v = 0.00001. k = y/70/3. 



Since L is negative the latter term will slow the accel- 
eration in the A — >• 0, A —> oo limit. Thus we conclude 
that for models with fixed H those with a larger value 
of A will provide a longer period of inflation. We want 
to stress that this conclusion is completely dependent on 
the holographic description through evolution associated 
with the p direction - we are learning about the role of the 
strong interaction dynamics. One certainly can not, for 
example, deduce the motion from the simplistic approx- 
imate potentials we displayed in Fig 3. More generally 
we note that the holographic dependence on p introduces 
considerably more complication to the evolution and ex- 
plicit simulation is required. 

It is clear that if we wish to prolong the early time roll 
period we need to reduce the rate of change of the dilaton 
and push A far above the vacuum energy H . There are 
two ways we can do that within our ansatz (21) - we can 
reduce A for a fixed V so that the step is smaller. We 
would expect this to increase the roll time. Also for fixed 
A we can decrease T so the change in A occurs over a 
larger p range. This reduces the dilaton derivative but 
spreads the region of p over which the change is occur- 
ring so there could be no net change to the total rate of 
acceleration - we will need to test this case numerically. 

In Fig 9 we plot e for the three different sets of param- 
eters in /?. Like Fig 5, e + 3p is always negative before 
it suddenly vanishes. Even though our computation is 
not valid when e starts changing fast we choose to take 
the time when s + ip vanish as the end time of inflation, 
t e . Since our purpose is to compare the inflation time for 
different dilaton profiles qualitatively, this approximation 
will not make a difference. 

We will now make a comparison of the roll time for a 
number of different dilaton step profiles. In each case the 
difference in vacuum energy between the symmetric and 
vacuum symmetry breaking configuration is the same. 
We also use the same initial velocity perturbation for 
the D7 brane (30). As an example of the differences we 
plot the energy density e against time in Fig 9 and the 
condensate c against time for three configurations in Fig 
10. First compare the curves for the dilaton parameters 
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FIG. 11: t e against A for a sequence of values of Fs: 
1.5(Black), l(Blue), 0.5 (Red), v = 0.00001. H = ^70/3. 



[A = 3, r = 1, A = 3.24] and [A = 1.8, T = 1, X = 4.325]. 
Decreasing A indeed increases the time the configuration 
takes to reach the tipping point to the true vacuum. Next 
we can change T to try to further increase the roll time - 
the final curve is for the configuration [A = 1.8, T = 0.5, 
A = 5.882] and indeed we find a further lengthening of 
the inflationary period. 

Finally we show these trends in more detail in Fig 11 
where wc plot t e against A for a sequence of values of T 
(1.5, 1, 0.5). There is a clear trend for decreasing both A 
and r increasing the roll time. As this roll time increases 
the parameters become more fine tuned reflecting the fine 
tuning wc arc making in the effective potential for c, ie 
the usual fine tuning in inflation. 

This analysis has been performed again in a heav- 
ily damped scenario which as discussed in (36) assumes 
a growing value of N in the gauge theory as the en- 
ergy density is reduced relative to the Planck scale. Wc 
choose this regime primarily for computational conve- 
nience. Having high damping removes late time oscil- 
lations of the D7 motion. It also extends the roll time 
making changes in that roll time more easily apparent. 
To demonstrate that the effects we have observed are still 
present at lower values of the damping we finally display 
the behaviour of the condensate against time for the case 
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FIG. 12: t e against A for two cases sequence of values. Blue: 
[A = 27.3, r = 0.5, A = 0.25], Red: [A = 2, T = 0.5, A = 
5.326]. v = 0.00001. H = y/l/3. 

H = y/l/Z and for steps with two different values of step 
height A in Fig 12. The late time oscillations are now ap- 
parent but the increase of the roll time with decreasing 
A is also maintained. 



V. DISCUSSION 

Inspired by small field inflation models, we have used 
holography to study out of equilibrium dynamics in 
a strongly coupled gauge theory with chiral symmetry 
breaking. We have shown that we can extract an ap- 
proximate effective potential for the quark condensate. 
We have then been able to explicitly follow the roll of 
a configuration from the symmetric potential maximum 
to the symmetry breaking potential minimum. We first 
performed this computation in a model with chiral sym- 
metry breaking induced by a magnetic field. While the 
configuration lies near the potential maximum it gener- 
ates inflation. Formally we worked in a probe limit that 
requires that the rate of inflation is small relative to the 
typical relaxation time of the background N = 4 gauge 
dynamics. 

We then studied a holographic model of a large iV c 
gauge theory with an arbitrary running of the gauge cou- 
pling. It is a model in the spirit of AdS / QCD since we do 
not backreact the geometry to the dilaton profile. The 
dilaton profile we chose to study means the gauge the- 
ory has a strongly coupled conformal UV and deforms to 
an IR conformal theory with a higher value of the gauge 
coupling. That change in the coupling triggers chiral 
symmetry breaking for the quarks. We have found that 
when the change in the coupling is fine tuned to be small 
or slow running the theory can give rise to inflationary 
dynamics. An important contribution to this slowing of 
the roll was made by the D7 dynamics in the radial di- 
rection of the geometry through (39) - this suggests that 
the strong dynamics of these gauge theories encourages 
inflation. Are there broad lessons to be learnt here about 
strongly coupled gauge theories or is this an oddity of a 
particular AdS set up? 



There is a growing belief that many asymptotically 
free gauge theories indeed give rise to strongly coupled 
conformal regimes. Scibcrg's dualities for 3Sf = 1 super- 
symmetric QCD [29] were the first hint - they show that 
SQCD flows to a non-trivial IR conformal theory in the 
range N c + 1 < Nf < 3N C . For much of that range 
the UV degrees of freedom are strongly coupled. Near 
Nf = 3N C and at large N c these phases match onto 
the perturbative Banks Zaks fixed points [30]. Banks 
Zaks fixed points also exist in non-supersymmetric the- 
ories and it is reasonable to expect that IR conformal 
fixed points also exist for a considerable range of Nf in 
those theories (see [31], [32] for some speculation about 
these theories). Theories with higher dimension repre- 
sentation matter fields would also be expected to gener- 
ate strongly coupled IR fixed points. Recently there has 
begun to be lattice simulations of QCD with varying un- 
quenched quark flavours [33, 34] and higher dimensional 
representation matter present [35, 36] - there is certainly 
encouragement in these results for the view that IR con- 
formal theories exist in some of these cases. 

If we believe that such fixed points are fairly common 
then we can imagine several ways for how to construct 
a theory with the profile for the running coupling we 
have studied in this paper. For example, one could begin 
with an SU(iV c ) theory with sufficient non-fundamental 
matter to place it at a strongly coupled fixed point. If 
N c is appropriately large then we can also add funda- 
mental quark multiplcts as a perturbation - their contri- 
butions to the beta function coefficients will be Nf/N c 
suppressed and hence they will most likely generate just 
a small change in the fixed point's coupling value. The 
fundamental quarks though can be used to dial a profile 
for the running coupling if they can be sequentially de- 
coupled. For example, if they are vector-like one could 
just put in masses to adjust the running at the order 
Nf/N c shifting the IR theory from one conformal point 
to another. This realizes the running we described above. 
It is not obvious how chiral symmetry breaking might 
happen in this scenario though. Usually one imagines 
an NJL-model-like critical coupling for chiral symmetry 
breaking [27]. One would need to tune the shift in the 
coupling to cross that critical value for the fundamental 
quarks. Naively higher dimension representation quarks 
might also be expected to condense at the same point or 
at a lower value of the coupling though. 

Another possibility, which is sometimes discussed [37], 
is that pure glue gauge theories in fact have an IR fixed 
point for the coupling. One could imagine perturb- 
ing that fixed point with some fundamental matter and 
again, by appropriate decoupling, change the running to 
match that we seek. The pure glue fixed point would 
again need to lie just above the critical coupling for chi- 
ral symmetry breaking. 

Finally the other possibility we are led to is walking dy- 
namics [38, 39]. A theory with fundamental quarks that 
is approaching an IR fixed point where the coupling is 
tuned just above the critical coupling might spend many 
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RG decades at strong coupling but without triggering chi- 
ral symmetry breaking before finally reaching the critical 
value at a low scale. Equally the coupling might cross the 
critical value at a high scale but be so close to the critical 
value from above that the quarks' dynamical mass is at 
a much lower scale so they don't decouple and the fixed 
point is maintained. Such a theory would naturally re- 
alize both the small A and small T limit of our coupling 
ansatz, both of which led towards inflation. Walking is 
of course proposed as a solution of the mini-hierarchy 
problem in technicolour theories of electroweak symme- 
try breaking [40, 41]. It is certainly intriguing if the fine 
tuning already used to solve that problem also generates 
inflationary dynamics. 

In conclusion it certainly seems possible that a range 



of asymptotically free gauge theories might realize the 
behaviour we have seen and generate effective dynamics 
that encourages inflation. The need for intrinsic fine tun- 
ing, as usual in small field inflation, seems unavoidable 
though. 
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